Method of obtaining a transmission gain function

ABSTRACT

Method of obtaining a transmission gain function by means of an array of antennae, a signal to be transmitted by the array being weighted by a vector ({overscore (b)} d ) of N complex coefficients, referred to as the transmission weighting vector, N being the number of antennae in the array, the array transmitting to a telecommunication terminal over a transmission channel, referred to as the downlink channel, a downlink transmission signal (S d ) and the said terminal transmitting to the said array over a transmission channel, referred to as the uplink channel, an uplink transmission signal (S u ), the said downlink channel being disturbed by an isotropic noise (N′) and/or a directional noise, referred to as the downlink interference (I d ), the said transmission weighting vector ({overscore (b)} d ) being determined by means of a matrix product from a noise power matrix (D d ) which is a function of the power of the said isotropic noise and/or of the power of the said directional noise and a vector ({overscore (C d )}), referred to as the downlink channel vector, representing an angular sampling of the transfer function of the downlink channel in M directions k, k=0, . . . , M−1, belonging to the angular range covered by the array.

[0001] The present invention concerns in general terms a method of obtaining a transmission gain function. More particularly, the present invention relates to a method of obtaining an antenna gain in transmission mode for a base station in a mobile telecommunication system. It makes it possible to obtain an antenna gain function in transmission mode from an antenna gain function in reception mode.

[0002] The formation of channels or the elimination of interfering signals is well known in the field of narrow-band antenna processing. Each uses an array of antennae, generally linear and uniform (that is to say with a constant pitch) and a signal weighting module. More precisely, if it is wished to form a channel in reception mode, the signals received by the different antennae are weighted by a set of complex coefficients before being added. Conversely, if it is wished to form a channel in transmission mode, the signal to be transmitted is weighted by a set of complex coefficients and the signals thus obtained are transmitted by the different antennae.

[0003]FIG. 1 illustrates a known device for obtaining an antenna gain in transmission and reception. The device comprises an array of antennae (10 ₀),(10 ₁), . . . , (10 _(N−1)), a transmission weighting module (11) and a reception weighting module (15). The signals received by the different antennae, (x_(i)), i=0 . . . N−1, are weighted at (13 ₀),(13 ₁), . . . , (13 _(N) ⁻¹) by a set of complex coefficients (b_(ui)), i=0, . . . , N−1 before being added at (14) in order to give a signal R_(u). Conversely, a signal to be transmitted S_(d) is weighted as (12 ₀), (12 ₁) . . . , (12 _(N−1)) by a set of complex coefficients (b_(di)), i=0, . . . , N−1, before being transmitted by the different antennae.

[0004] If the vector of the received signals and the vector of the weighting coefficients are respectively denoted {overscore (x)}=(x₀, x₁ . . . , x_(N−1))^(T) and {overscore (b)}_(u)=(b_(u0), b_(u1) . . . , b_(uN−1))^(T), it is possible to write: $\begin{matrix} {R_{u} = {{\overset{\_}{b}}_{u}^{T}\overset{\_}{x}}} & (1) \end{matrix}$

[0005] The complex gain (or the complex gain function of the antenna) in reception mode can be written: $\begin{matrix} {{G\left( {\overset{\_}{b_{u}},\theta} \right)} = {{{\overset{\_}{b_{u}}}^{T} \cdot \overset{\_}{e_{u\quad \theta}}} = {\sum\limits_{i = 0}^{N - 1}\quad {b_{ui} \cdot {\exp \left( {{- {j\phi}}\quad i} \right)}}}}} & (2) \end{matrix}$

[0006] where {overscore (e_(uθ))} represents the vector {overscore (x)} corresponding to a flat wave arriving at an angle of incidence θ, and

φ_(i)=(2πd/λ).i. sin(θ)=(2πdf/c).i. sin(θ)  (3)

[0007] is the difference in operation between consecutive antennae for a uniform linear array with a pitch d, λ and f being respectively the wavelength and the frequency of the flat wave in question; and

φ_(i)=2πRΔθ/λ. sin(θ−θ_(i))=2πRfΔθ/c. sin(θ−θ_(i))  (4)

[0008] for a circular array where θ_(i) is the angle between a reference axis and the normal to the antenna of index i, R the radius of curvature of the array, Δθ is the angular difference between two consecutive antennae in the array.

[0009] Likewise, the complex gain (or the complex gain function) in transmission can be written: $\begin{matrix} {{G\left( {\overset{\_}{b_{d}},\theta} \right)} = {{{\overset{\_}{b_{d}}}^{T} \cdot \overset{\_}{e_{d\quad \theta}}} = {\sum\limits_{i = 0}^{N - 1}\quad {b_{di} \cdot {\exp \left( {{j\phi}\quad i} \right)}}}}} & (5) \end{matrix}$

[0010] with the same conventions as those adopted above and where {overscore (e_(dθ))} represents the vector {overscore (x)} corresponding to a flat wave transmitted in the direction θ.

[0011] The weighting vectors in reception and transmission mode will be called respectively {overscore (b_(u))} and {overscore (b_(d))}.

[0012] When the array of antennae is functioning in reception mode at a given frequency, different known methods, notably the Wiener filtering method, make it possible to determine the weighting vector {overscore (b_(u))} which maximises the signal to noise ratio. In a mobile telecommunications system, the array of antennae of a base station receives signals transmitted by a plurality of mobile terminals. In the context of a transmission in CDMA (Code Division Multiple Access) mode, the signals transmitted by the different mobile terminals are separated by means of the use of orthogonal codes on transmission and filters adapted to these codes on reception. In practice, however, the separation of the different signals received is not perfect. For an uplink between a given mobile terminal and the base station which serves it, the criterion to be maximised is then the ratio of signal to noise plus interference, the latter being due to the signals transmitted by the other mobile terminals. Likewise, the downlink between a base station and a given mobile terminal is disturbed not only by the background noise but by the interference due to the signals transmitted by the said base station to other mobile terminals. Though it is relatively easy to optimise the weighting vector in reception mode, {overscore (b_(u))}, by estimating the uplink channel and the density of interference at the base station, it is quite different with regard to the optimisation of the weighting vector in transmission mode, {overscore (b_(d))}. This is because the estimation of the downlink channel and the density of interference cannot be made directly at the base station and a transmission of this information by the mobile terminals is necessary. However, this transmission of information consumes conveyance resources on the uplink, which can be disadvantageous, notably in the case of rapid variations in the channel transfer function, for example when the mobile terminal is moving at high speed.

[0013] The aim of the invention is to propose a method of determining the transmission weighting vector, {overscore (b_(d))}, optimising the ratio of signal to noise plus interference on the downlink and requiring the transmission only of a small quantity of information on the uplinks.

[0014] To this end, the invention is defined by a method of obtaining a gain function on transmission by means of an array of antennae, a signal to be transmitted by the network being weighted by a vector ({overscore (b)}_(d)) of N complex coefficients, referred to as a transmission weighting vector, N being the number of antennae in the array, the array transmitting to a telecommunication terminal on a transmission channel, referred to as the downlink, a downlink transmission signal (S_(d)), and the said terminal transmitting to the said network on a transmission channel, referred to as the uplink channel, an uplink transmission signal (S_(u)), the said downlink channel being interfered with by an isotropic noise (N′) and/or a directional noise, referred to as the downlink interference (I_(d)), the said transmission weighting vector ({overscore (b)}_(d)) being determined by means of a matrix product from a noise power matrix (D_(d)) which is a function of the power of the said isotropic noise and/or the power of the said directional noise and a vector ({overscore (C_(d))}), the said downlink channel vector, representing an angular sampling of the transfer function of the downlink channel in M directions k, k=0, . . . , M−1, belonging to the angular range covered by the array.

[0015] Advantageously, the said downlink channel vector ({overscore (C_(d))}) is obtained from variations in the transfer function of the uplink channel. The said downlink channel vector ({overscore (C_(d))}) is obtained for example from variations ({overscore (ΔC_(u))}) in a vector ({overscore (C_(u))}), referred to as the uplink channel vector, representing an angular sampling of the transfer function of the uplink channel in the said M directions.

[0016] The variations ({overscore (ΔC_(d))}) in the downlink channel vector can be obtained from variations ({overscore (ΔC_(u))}) in the uplink channel.

[0017] The variations Δc_(dk) in the components c_(dk) of the downlink channel vector ({overscore (C_(d))}) are advantageously obtained by means of variations Δc_(uk) in the components c_(uk) of the uplink vector by: Δc_(dk)/c_(dk)=f_(d)/f_(u).Δc_(uk)/c_(uk) whereby f_(u) is the frequency used on the said uplink channel and f_(d) is the frequency used on the said downlink channel.

[0018] The said downlink channel vector ({overscore (C_(d))}) is then obtained by integration of the said variations ({overscore (ΔC_(d))}) in the said downlink channel vector and an initial value ({overscore (C_(d))}(0)) transmitted by the said terminal.

[0019] According to a first embodiment, the noise matrix is a diagonal matrix of size M×M and components {square root}{square root over (σ_(dk) ²+γ_(d)N′₀/I_(d))} where σ_(dk) ² is the downlink interference power in the direction k, N′₀ is the power of the isotropic noise, γ_(d)=1/∥{overscore (C)}_(d)∥² and I_(d) is the total power of the downlink interference.

[0020] According to a second embodiment, the array transmitting on a plurality of downlink channels a plurality of transmission signals to a plurality of telecommunication terminals and receiving from them a plurality of transmission signals transmitted on a plurality of uplink channels, each downlink channel j relating to a terminal j of the said plurality being associated with a transmission weighting vector {overscore (b)}_(d)(j), the second noise matrix relating to the downlink channel j is a diagonal matrix of size M×M and components {square root}{square root over (σ_(dk) ²(j)+γ_(d)(j).N′₀/I_(d)(j))} where σ_(dk) ²(j) is the power of the downlink interference for the downlink channel j in the direction k, γ_(d)(j) is a coefficient characterising the power transfer on the downlink channel j, N′₀ is the power of the second isotropic noise, and I_(d) is the total power of the downlink interference.

[0021] Advantageously, the coefficient γ_(d) (j) is transmitted to the array by the terminal j on the associated uplink channel.

[0022] Alternatively, the coefficient γ_(d) (j) is estimated by the base station from a coefficient (Γ) characterising the power transfer in the uplink direction.

[0023] For a given downlink channel j, the downlink interference power in the direction k, σ_(dk) ²(j), can be estimated according to the power of the signals (S_(d)(j′)) transmitted on the downlink channels j′ distinct from j, a coefficient β_(d)(j) characterising the orthogonality of the downlink channel j, components g_(dk)(j′)) of the gain vectors ({overscore (G)}_(d)(j′)) relating to the said distinct downlink channels j′, the gain vectors consisting of an angular sampling in the said M directions of the transmission gain functions obtained for the said distinct downlink channels j′.

[0024] Finally, the said coefficient β_(d)(j) can be estimated from a coefficient characterising the orthogonality of the uplink channel j.

[0025] The invention is also defined by a device adapted to implement the method disclosed above.

[0026] The characteristics of the invention mentioned above, as well as others, will emerge more clearly from a reading of the following description given in relation to the accompanying figures, amongst which:

[0027]FIG. 1 depicts schematically a known device for obtaining an antenna gain function;

[0028]FIG. 2 depicts schematically an uplink transmission channel in a mobile telecommunication system;

[0029]FIG. 3 depicts schematically a downlink transmission channel in a mobile telecommunication system;

[0030]FIG. 4 depicts schematically a device for obtaining an antenna gain function according to one embodiment of the invention.

[0031] A first general idea at the basis of the invention is to sample the transmission and reception gain functions in order to construct transmission and reception gain vectors. As will be shown, optimum weighting vectors, in terms of ratio of signal to noise plus interference, can then be obtained from transmission and reception gain vectors according to matrix equations.

[0032] A second general idea at the basis of the invention is to obtain a transmission weighting vector, optimum in terms of ratio of signal to noise plus interference obtained, according to the reception gain weighting vector assumed itself to be optimum.

[0033] A third general idea at the basis of the invention is to estimate the downlink channel from variations in the uplink channel.

[0034] It will shown first of all that the weighting vectors can be obtained from a series of samples of the corresponding gain function.

[0035] Consider first of all a linear uniform array, formed by N antennae spaced apart by a pitch d and operating at the frequency f. The antennae gain function G₀(θ), obtained in the absence of any weighting (that is to say with {overscore (b₀)}=(1, 1, . . . , 1)^(T)): $\begin{matrix} {{{G_{0}(\theta)}} = {{\frac{\sin \left( {N\quad {\phi/2}} \right)}{\sin\left( \quad {\phi/2} \right)}\quad {with}\quad \phi} = {2\pi \quad {{{fd}/c} \cdot \sin}\quad \theta}}} & (6) \end{matrix}$

[0036] This function has zeros for the values σ_(k)=2kπ/N, k a non-zero integer such that σ_(k)ε[−π,π[, that is to say in the directions for which sinθ_(k)=k.c/Nfd, when this expression has a direction. The phase difference between two consecutive zeros in the gain diagram is constant and is equal to Δσ=2π/N. The angular difference between two consecutive zeros of the diagram varies in terms of Arc sin., a function whose derivative is increasing on [−1,1] and is therefore minimum for the angular difference between the first and second zeros. It is therefore bounded by Δθ_(min)=c/Nfd if N is sufficiently great. It will be assumed that the frequencies used are less than f₀ where f₀ is the natural frequency of the array. It can be concluded from this that the spectrum of the function G₀(θ) has a support bounded by 1/Δθ_(min)=N/2.

[0037] In more general terms, let G(θ) be the antenna gain function obtained by means of a weighting vector {overscore (b)}. G can be expressed as the Fourier transform (FT) (in reception) or the inverse Fourier transform (in transmission) of the complex weighting distribution of the antenna, namely: ${b(x)} = {\sum\limits_{i = 0}^{N - 1}\quad {b_{i} \cdot {\delta \left( {x - x_{i}} \right)}}}$

[0038] with x_(i)=i.d; this gives:

[0039] G_(u)(θ)=B(sinθ) with B(u) = ∫_(−∞)^(+∞)b(x)exp (−j2πux/λ)  x

[0040] and likewise G_(d)(θ)=B′(sinθ) with B^(′)(u) = ∫_(−∞)^(+∞)b(x)exp (j2πux/λ)  x.

[0041] The function b(x) being bounded by N.d, the difference between two zeros of the function B or B′ is at least λ/N.d and therefore all the more so 2/N. Given the increase in the derivative of the function Arc sin. the minimum difference between two zeros of the function G is 2/N. The function G therefore has a spectrum bounded by N/2.

[0042] According to the Shannon sampling theorem, it is concluded from this that it is possible to reconstitute the function G(θ) if sampling is carried out at a frequency greater than the Nyquist frequency, that is to say N. In other words, for an angular range [−π/2,π/2], a minimum of M>π.N samples are necessary, where M is an integer. In practice K.N samples can be taken with K integer, K≧4.

[0043] For a circular array, it can be shown that 1/Δθ_(min)=N and the angular range being [−π,π], M (M>π.N and M integer) angularly equidistributed samples also suffice to reconstitute the function G(θ).

[0044] In the general case of the sampling of any gain function G(θ), it is necessary first to filter G(θ) by means of an anti-aliasing filter before sampling it. It then suffices to take M samples of the diagram filtered over the entire angular range in order to reconstitute the filtered diagram.

[0045] The samples of a gain function, possibly filtered by anti-aliasing filtering, if necessary, will be denoted g_(k), k=0, . . . , M−1 that is to say g_(k)=G′(θ_(k)) where the instances of θ_(k) are M angles equidistributed over [−π/2,π/2] or [−π,π] and where it is assumed that G′ was the filtered version of the reference complex diagram.

[0046] It is now possible to define a linear application, h^(f) _(s) of C^(N) in C^(M) which makes the vector h_(s) ^(f)({overscore (b)})={overscore (G)}=(g₀, g₁, . . . , g_(M−1))^(T) where g_(k)=G({overscore (b)},θ_(k)) correspond to any weighting vector {overscore (b)}. The image of C^(N) by h^(f) _(s), is a vector sub-space of C^(M) with a dimension at most equal to N, which will be denoted Im_(f). If a base of C^(N) is chosen, for example the canonical base, and a base of C^(M), it is possible to express the linear application h^(f) _(s) by a matrix H_(f) of size M×N which is at most of rank N.

[0047] Let {overscore (G)} be any gain vector corresponding to a sampled gain function. Seek the vector {overscore (b)} such that h_(s) ^(f)({overscore (b)}) is as close as possible to {overscore (G)} in the sense of a certain metric. There will be taken as a norm the Euclidian norm on C^(M), namely ${\overset{\_}{G}}^{2} = {\sum\limits_{k = 0}^{M - 1}\quad {{g_{k}}^{2}.}}$

[0048] If it exists, the sought-for vector {overscore (b)} is then such that h_(s) ^(f)({overscore (b)})={overscore (G)}_(p) where {overscore (G)}_(p) is the orthogonal projection of the vector {overscore (G)} onto Im_(f). If the matrix H_(f) is of rank N, the sought-for vector {overscore (b)} exists and can be written:

{overscore (b)}=H _(f) ⁺ .{overscore (G)}  (7)

[0049] where H_(f) ⁺=(H_(f)*^(T).H_(f))⁻¹.H_(f)*^(T) is the pseudo-inverse matrix of the matrix H_(f) with H_(f)*^(T) a conjugate transpose of the matrix H_(f).

[0050] In order to express the matrix H_(f), it is necessary to agree on a base of the starting space and a base of the arrival space. It is possible to choose as the base of C^(M) the canonical base and as a base of C^(N) a base adapted to the description of the flat waves of frequency f. Consider the distinct vectors {overscore (e)}_(k), k=0, . . . , N−1, such that {overscore (e)}_(k)=(e_(k,0), e_(k,1), . . . , e_(k,N−1))^(T) with $e_{k,i} = {{\exp \left( {{j \cdot \frac{2\pi \quad {fd}}{c} \cdot  \cdot \sin}\quad \theta_{k}} \right)} = {\exp \left( {j\quad {\pi \cdot \eta \cdot  \cdot \sin}\quad \theta_{k}} \right)}}$

[0051] with η=f/f₀ and where the instances of θ_(k) belong to the interval [−π/2, π/2]. The vectors {overscore (e)}_(k) are the weighting vectors of the array making it possible to form beams in the directions θ_(k). The vectors {overscore (e)}_(k) form a base if the determinant of the coordinates of the instances of {overscore (e)}_(k) in the canonical base of C^(N) is non-zero. This determinant is a Vandermonde determinant which is equal to $\prod\limits_{p \neq q}\left( {{\exp \left( {j\quad \phi_{p}} \right)} - {\exp \left( {j\quad \phi_{q}} \right)}} \right)$

[0052] with φ_(k)=πη sin θ_(k). This determinant is cancelled out if and only if there are two angles θ_(p) and θ_(q) such that sin θ_(p)−sin θ_(q)=2/η. In other words, for η<1 the N vectors {overscore (e)}_(k) always form a base and for η=1 only the case θ_(p)=−θ_(q)=π/2 is excluded. The directions can, for example, be chosen so as to be equally distributed, that is to say such that θ_(k)=kπ/N with k=−(N−1)/2, . . . , 0, . . . , (N−1)/2. In this case, the matrix H_(f) has as its components: $H_{pq} = {\sum\limits_{i = 0}^{N - 1}\quad {{\exp \left( {{j\pi}\quad {\eta \cdot  \cdot {\sin \left( {p\quad {\pi/N}} \right)}}} \right)}{\exp \left( {{- j}\quad {{\pi\eta} \cdot  \cdot {\sin \left( {q\quad {\pi/M}} \right)}}} \right)}}}$

[0053] or: $\begin{matrix} {H_{pq} = {{\sum\limits_{i = 0}^{N - 1}\quad {\exp \left( {j\quad \pi \quad {\eta \cdot  \cdot \left\lbrack {{\sin \left( {p\quad {\pi/N}} \right)} - {\sin \left( {q\quad {\pi/M}} \right)}} \right\rbrack}} \right)}} = {{\exp \left( {{j\left( {N - 1} \right)}{\Psi_{pq}/2}} \right)} \cdot \frac{\sin \left( {N\quad {\Psi_{pq}/2}} \right)}{\sin \left( {\Psi_{pq}/2} \right)}}}} & (8) \end{matrix}$

[0054] with Ψ_(pq)=πη(sin(pπ/N)−sin(qπ/M)).

[0055] Alternatively, it is possible to choose as a starting base another base adapted to the frequency f, the one formed by the vectors {overscore (e)}′_(k), such that e′_(′k,i)=exp(jπ.η.i. sinθ_(k)) with sinθ_(k)=2k/ηN and k=−(N−1)/2, . . . , 0, . . . , (N−1)/2. These vectors exist if |sin θ_(k)|≦1,∀k, that is to say for η>1−1/N, and in this case the vectors {overscore (e)}′_(k) form a base which has the advantage of being orthogonal.

[0056] Alternatively, it is possible to choose as the starting base the canonical base of C^(N), which has the advantage of not depending on the frequency. In this case, the matrix H′_(f) expressed in this base is written:

H′ _(f) =H _(f) .T ⁻¹  (9)

[0057] where T is the matrix of the coordinates of {overscore (e)}_(k) in the canonical base, that is to say T_(pp′)=exp(jπp sin(p′/N)). It was seen above that this matrix had a non-zero Vandermonde determinant and was consequently not reversible.

[0058] Assume now that it is sought to approximate a gain function obtained at a first frequency f₁, f₁ ≤ f₀

[0059] and denote the vector of the samples associated with this gain function.

[0060] Let a second working frequency be f₂, f₂≦f₀. {overscore (G₁)} belonging to C^(M), if the matrix H_(f2) is of rank N, it is possible to find a vector {overscore (b₂)} such that h_(s) ^(f2)({overscore (b₂)}) is the projection of h_(s) ^(f1)({overscore (b₁)}) onto Imf₂. The vector {overscore (b₂)} is obtained by means of the matrix equation:

{overscore (b₂)}= H _(f2) ⁺ .H _(f1) {overscore (b₁)}  (10)

[0061] This equation makes it possible in particular to obtain, at a second working frequency, a sampled gain diagram which is as close as possible to the one, referred to as the reference, obtained at a first working frequency.

[0062] Equation (10) applies advantageously to the array of antennae of a base station in a mobile telecommunication system operating in FDD (Frequency Division Duplex) mode. In such a system, a frequency f_(d) is used on the downlinks and a frequency f_(u) distinct from f_(d) is used on the uplinks. Equation (10) then makes it possible to directly obtain the weighting vector on transmission {overscore (b_(d))} from the weighting vector on reception {overscore (b_(u))}.

{overscore (b_(d))}= H _(d) ⁺ .H _(u) {overscore (b_(u))}  (11)

[0063] where there has been noted: H_(d)=H_(fd) and H_(u)=H_(fu).

[0064] Equation (11) makes it possible, as has been seen, to obtain, at the transmission frequency f_(u), a sampled gain diagram which is as close as possible to a reference diagram obtained at the reception frequency f_(d). However, the interference profile, that is to say the angular distribution of the power of the interference, is not necessarily the same on the downlink channel as on the uplink channel. This is because the directions of the interfering sources are not necessarily identical in transmission and reception. Consequently, though the reception gain diagram is optimum for a reception interference profile, it will not necessarily be so for a transmission interference profile. As will be shown later, if the transmission and reception interference profiles differ, equation (11) must be modified in order to take account of this difference.

[0065]FIG. 2 depicts the assembly consisting of the uplink channel (20), the array of antennae (22) and the reception weighting module (23). The effect of the noise has been represented by the addition (21) of a directional noise {overscore (I)}_(u) due to the interfering signals, and at (24) an isotropic centred white Gaussian background noise N.

[0066] Just as the gain function can be represented by a gain vector, the channel can be modelled by a dimension vector M, defined as the angular sampling of the transfer function of the channel in the directions θ_(k), k=0, . . . , M−1 and denoted {overscore (C)}_(u)=(c_(u0), c_(u1), . . . , c_(uM−1))^(T). This vector has P amongst M non-zero coefficients, where P is the number of channel propagation paths. For these P coefficients c_(uk), c_(uk)−α_(uk).exp−j(2πf_(u).L_(uk)/c+φ_(uk)) where L_(uk) is the length of the path concerned, α_(uk) the coefficient of attenuation of the signal propagating along the said path and φ_(uk) the polarisation of the incident signal.

[0067] The signal R_(u) received by the base station can be written: $\begin{matrix} {R_{u} = {{{\overset{\_}{G}}_{u}^{T} \cdot \left( {{S_{u} \cdot {\overset{\_}{C}}_{u}} + {\overset{\_}{I}}_{u}} \right)} + N}} & (12) \end{matrix}$

[0068] where {overscore (G)}_(u) is the reception gain vector and S_(u) is the signal transmitted by the mobile terminal.

[0069] The ratio of signal to noise plus interference is: $\begin{matrix} {\left( {{C/I} + N} \right)_{u} = {\frac{E\left( {{{{\overset{\_}{G}}_{u}^{T} \cdot S_{u}}{\overset{\_}{C}}_{u}}}^{2} \right)}{{E\left( {{{\overset{\_}{G}}_{u}^{T} \cdot {\overset{\_}{I}}_{u}}}^{2} \right)} + {E\left( {N}^{2} \right)}} = \frac{P_{u} \cdot {{{\overset{\_}{G}}_{u}^{T} \cdot {\overset{\_}{C}}_{u}}}^{2}}{N_{0} + {I_{u} \cdot {\sum\limits_{k = 0}^{M - 1}\quad {\sigma_{uk}^{2}{g_{uk}}^{2}}}}}}} & (13) \end{matrix}$

[0070] where P_(u) is the power of the signal S_(u), N₀ is the power of the background noise and {overscore (I)}_(u)={square root}{square root over (I_(u))}.{overscore (I)}_(u) ⁰ where {overscore (I)}_(u) ⁰ is the standardised vector whose components are assimilated to standardised centred Gaussian random variables N(0,σ_(uk) ²) that is to say such that ${\sum\limits_{k = 0}^{M - 1}\quad \sigma_{uk}^{2}} = 1$

[0071] and where I_(u)=E(∥{overscore (I)}_(u)∥²) is the total power of the directional noise (that is to say of the interference on the uplink channel). It was assumed in (13) that the isotropic noise was independent of the directional noise.

[0072] Expression (13) can then be written: $\begin{matrix} {\left( {{C/I} + N} \right)_{u} = {\frac{P_{u}}{I_{u}} \cdot \frac{{{{\overset{\_}{\Lambda}}_{u}^{T}{\overset{\_}{\Omega}}_{u}}}^{2}}{{{\overset{\_}{\Lambda}}_{u}}^{2}}}} & (14) \end{matrix}$

[0073] where Ω_(u)=D_(u) ⁻¹.C_(u), Λ_(u)=D_(u)G_(u) and D_(u)=Diag({square root}{square root over (σ_(uk) ²+γ_(u)N₀/I_(i))}) with γ_(u)=1/∥{overscore (G)}_(u)∥². Expression (14) is maximal for {overscore (Λ)}_(u)=Ω _(u)* and therefore for:

{overscore (G)} _(u) =D _(u) ⁻¹{overscore (Ω)}_(u) *=D _(u) ⁻² .{overscore (C)} _(u)*  (15)

[0074] The reception weighting vector, {overscore (G_(u))}, optimum in the sense of the maximisation of the ratio of signal to noise plus interference on the channel, can then be expressed:

{overscore (b)} _(u) =H _(u) ⁺ .D _(u) ⁻¹{overscore (Ω)}_(u) *=H _(u) ⁺ .D _(u) ⁻² {overscore (C)} _(u)*  (16)

[0075]FIG. 3 depicts the assembly consisting of the downlink channel (30), the array of antennae (32) and the transmission weighting module (33). The effect of the noise has been represented by the addition at (31) of a directional noise {overscore (I)}_(d) due to the interfering signals and at (34) by a centred isotropic white Gaussian background noise N′.

[0076] Just like the uplink channel, the downlink channel can be modelled by a vector of dimension M, defined as the angular sampling of the transfer function of this channel in the directions θ_(k), k=0, . . . , M−1 and denoted {overscore (C)}_(d)=(c_(d0), c_(d1), . . . , c_(dM−1))^(T). This vector has P′ from amongst M non-zero coefficients, where P′ is the number of propagation paths of the channel. For these P′ coefficients c_(dk), c_(dk)=α_(dk).exp−j(2πf_(d).L_(dk)/c+φ_(dk)) where L_(dk′) is the length of the path concerned, α_(dk′) the coefficient of attenuation of the signal propagating on the said path and φ_(dk′) the polarisation of the incident signal.

[0077] The signal R_(d) received by the mobile terminal can be written: $\begin{matrix} {R_{d} = {{{\overset{\_}{C}}_{d}^{T} \cdot \left( {{S_{d} \cdot {\overset{\_}{G}}_{d}} + {\overset{\_}{I}}_{d}} \right)} + N^{\prime}}} & (17) \end{matrix}$

[0078] where {overscore (G)}_(d) is the transmission gain vector and S_(d) is the signal transmitted by the base station.

[0079] The ratio of signal to noise plus interference is: $\begin{matrix} {\left( {{C/I} + N} \right)_{d} = {\frac{E\left( {{{{\overset{\_}{C}}_{d}^{T} \cdot S_{d}}{\overset{\_}{G}}_{d}}}^{2} \right)}{{E\left( {{{\overset{\_}{C}}_{d}^{T} \cdot {\overset{\_}{I}}_{d}}}^{2} \right)} + {E\left( {N^{\prime}}^{2} \right)}} = \frac{P_{d} \cdot {{{\overset{\_}{C}}_{d}^{T} \cdot {\overset{\_}{G}}_{d}}}^{2}}{N_{0}^{\prime} + {I_{d} \cdot {\sum\limits_{k = 0}^{M - 1}\quad {\sigma_{dk}^{2}{{Cdk}}^{2}}}}}}} & (18) \end{matrix}$

[0080] where P_(d) is the power of the signal S_(d), and N′₀ is the power of the background noise, and where {overscore (I)}_(d)={square root}{square root over (I_(d))}.{overscore (I)}_(d) ⁰ where {overscore (I)}_(d) ⁰ is the standardised vector whose components are assimilated to standardised centred random Gaussian variables N(0,σ_(dk) ²), that is to say such that ${\sum\limits_{k = 0}^{M - 1}\quad \sigma_{dk}^{2}} = 1$

[0081] and where I_(d)=E(∥{overscore (I)}_(d)∥²) is the total power of the directional noise (that is to say of the interference on the downlink channel). It was assumed in (18) that the isotropic noise was independent of the directional noise.

[0082] Expression (18) can also be written: $\begin{matrix} {\left( {{C/I} + N} \right)_{d} = {\frac{P_{d}}{I_{d}} \cdot \frac{{{{\overset{\_}{\Lambda}}_{d}^{T}{\overset{\_}{\Omega}}_{d}}}^{2}}{{{\overset{\_}{\Lambda}}_{d}}^{2}}}} & (19) \end{matrix}$

[0083] where {overscore (Ω)}_(d)=D_(d) ⁻¹.{overscore (G)}_(d), Λ_(d)=D_(d){overscore (C)}_(d) and D_(d)=Diag({square root}{square root over (σ_(dk) ²+γ_(d)N′₀/I_(d))}) with γ_(d)=1/∥{overscore (C)}_(d)∥². Expression (19) is maximum for {overscore (Ω)}_(d)={overscore (Λ)}_(d)* and therefore for

{overscore (G)} _(d) =D _(d){overscore (Λ)}_(d) *=D _(d) ² {overscore (C)} _(d)*  (20)

[0084] The transmission weighting vector, {overscore (b_(d))}, optimum in the sense of the maximisation of the ratio of signal to noise plus interference on the downlink channel, can then be expressed as:

{overscore (b)} _(d) =H _(d) ⁺ .D _(d){overscore (Λ)}_(d) *=H _(d) ⁺ .D _(d) ² {overscore (C)} _(d)*  (21)

[0085] If it is assumed that the transfer function of the downlink channel is identical to that of the uplink channel, that is to say if {overscore (C_(d))}={overscore (C)}_(u), it is possible to derive from (16) and (21) the relationship between the optimum weighting vectors {overscore (b_(d))} and {overscore (b_(u))}:

{overscore (b)} _(d) =H _(d) ⁺ .D _(d) ² .D _(u) ² .H _(u) .{overscore (b)} _(utm ()21)

[0086] It should be noted that equation (11) is a particular case of equation (22) when D_(d)=D_(u) ⁻¹. This will be the case in particular if the noise on the uplink channel and the noise on the downlink channel consist solely of isotropic noise.

[0087] The matrix D_(u) can be estimated at the base station from a measurement of the noise power and the interference in the directions θ_(k), for example during a period of silence of the mobile terminal. On the other hand the matrix D_(d) cannot be estimated as simply.

[0088] Let it be stated that D_(d)=Diag({square root}{square root over (σ_(dk) ²+γ_(d)N′₀/I_(d))}) with

γ_(d)=1/∥{overscore (C)} _(d)∥²  (23)

[0089] γ_(d) can be estimated by the mobile terminal and transmitted over the uplink channel to the base station. Since γ_(d) changes only slowly over time, the quantity of information to be transmitted relating to this parameter will be small.

[0090] Advantageously, it will be assumed that the power transfer coefficient does not depend on the frequency and is identical for the downlink channel and the uplink channel, namely ∥{overscore (C_(d))}∥²=∥{overscore (C_(u))}∥²=1/Γ. This assumption will in particular be verified if {overscore (C_(d))}={overscore (C)}_(u). The value of Γ can then be estimated directly by the base station, for example at the power control loop.

[0091] It is also possible to estimate the power of the interference in the direction θ_(k), that is to say I_(d).σ_(dk) ², by expressing the fact that it is due to the transmission in the direction k of signals intended for mobile terminals TS_(j) other than the one considered, that is to say TS_(j0). Because of the overlap of the transmission gain diagrams for the signals intended for the mobile terminals TS_(j) on the one hand and the mobile terminal TS_(j0) on the other hand and the lack of orthogonality between these signals, the interference power attributable to the mobile terminals TS_(j) in the direction θ_(k) can be written: $\begin{matrix} {{{I_{d}\left( j_{0} \right)} \cdot {\sigma_{dk}^{2}\left( j_{0} \right)}} = {{\beta_{d}\left( j_{0} \right)} \cdot {\sum\limits_{j \neq j_{0}}{{{S_{d}(j)}}^{2} \cdot {{g_{dk}(j)}}^{2}}}}} & (24) \end{matrix}$

[0092] where the indices between parentheses have been added so as to distinguish the quantities relating to the different downlink channels (that is to say intended for the different mobile terminals) and where:

[0093] β_(d)(j₀) is the coefficient of orthogonality of the downlink channel, destined for TS_(j0);

[0094] S_(d)(j) is the power of the transmission signal destined for the terminal TS_(j);

[0095] g_(dk)(j) is the kth coefficient of the gain vector {overscore (G)}_(d)(j) relating to the transmission to TS_(j).

[0096] If it is assumed that the coefficient of orthogonality of the downlink channel, β_(d)(j₀), is little different from that of the uplink channel, β_(u)(j₀), the above three quantities are available at the base station without a return of information by the mobile terminal being necessary. As seen above, the power transfer coefficient, γ_(d)(j₀), is transmitted to the base station on the uplink channel from TS_(j0) or directly estimated by the latter. It is therefore possible to obtain the matrix D_(d) for a slight additional cost in terms of conveyance resources.

[0097] The only quantity in equation (24) liable to vary rapidly over time is the power of the transmission signals S_(d)(j). In the case of a transmission in DS-CDMA mode, it will for example be possible to update these power values at each transmission slot.

[0098] Equation (22) was obtained by assuming {overscore (C_(d))}={overscore (C)}_(u). This equation is not in general satisfied, notably because of the fact that the frequency used on the downlink channel is different from that used on the uplink channel. As a result the vector {overscore (b_(d))} calculated by means of equation (22) does not maximise the ratio of signal to noise plus interference, on reception by the mobile terminal.

[0099] The method according to the present invention proposes to calculate the optimum vector {overscore (b_(d))} from equation (21) and an estimation of the downlink channel, that is to say of the vector {overscore (C_(d))}. It should be stated that the vectors {overscore (C)}_(u)=(c_(u0), c_(u1), . . . , c_(uM−1))^(T) and {overscore (C)}_(d)=(c_(d0), c_(d1), . . . , c_(dM−1))^(T) consist of the transfer functions of the respectively uplink and downlink channels, sampled in the directions θ_(k), k=0, . . . , M−1. The vector {overscore (C_(u))} has non-zero components in the directions where the uplink channel has propagation paths. More precisely, if a propagation path exists in the direction θ_(k), c_(uk)=α_(uk).exp−j(2πf_(u).L_(uk)/c+φ_(uk)) where L_(uk) is the length of the path concerned, α_(uk) the coefficient of attenuation of the signal propagating along the said path, φ_(uk) is the polarisation of the incident signal and c_(uk)=0 otherwise. Likewise, the vector {overscore (C_(d))} has non-zero components in the directions where the downlink channel has propagation paths. More precisely, if a propagation path exists in the direction θ_(k′), c_(dk).α_(dk′).exp−j(2πf_(d).L_(dk)/c+φ_(dk)) where L_(dk′) is the length of the path concerned, α_(dk′) the coefficient of attenuation of the signal propagating along the said path, φ_(dk′) is the polarisation of the incident signal and c_(dk′)=0 otherwise. In the remainder of the explanations it will be assumed that the propagation paths of the uplink channel and of the downlink channel are identical, or in other words that L_(dk)=L_(uk).

[0100] The components of the vector {overscore (C_(u))}, c_(uk), can, for example, be determined, in a manner known per se, by means of pilot symbols transmitted by the mobile terminal. Advantageously, an estimation will be carried out of the coefficients of attenuation α_(uk) and of the directions of arrival of the paths conjointly, as described in the French patent application No. 00 11160 filed on 29.8.2000 in the name of the applicant. The vector {overscore (C_(u))}, relating to a given uplink channel, is first of all differentiated, or in other words the variation in the vector {overscore (C_(u))} during the interval of time Δt separating two consecutive estimations is evaluated. The variation in this vector can be written {overscore (ΔC)}_(u) =(Δc_(u0), Δc_(u1), . . . , Δc_(uM−1))^(T) where $\begin{matrix} {{\Delta \quad c_{uk}} \approx {{{- j} \cdot \left( {2\pi \quad f_{u}{L_{uk}/c}} \right) \cdot c_{uk} \cdot \frac{\partial L_{uk}}{\partial t}}\Delta \quad t}} & (25) \end{matrix}$

[0101] assuming that the coefficients of attenuation, α_(uk), of the different paths vary little in the course of the interval of time Δt. The vector {overscore (Δ)}c_(u) is next multiplied by the matrix M=Diag(f_(d)/f_(u).1/c_(uk)). The vector thus obtained, denoted {overscore (Δ)}_(d)=(δ₀, δ₁, . . . , δ_(M−1)) is then used for integrating the components of the vector {overscore (C_(d))} by virtue of the recurrence formula:

c _(dk)(t+Δt)=c _(dk)(t)(1+δ_(k)(t).Δt)  (26)

[0102] The initialisation of the calculation is effected by means of a vector {overscore (C_(d))}(0) of components c_(dk)(θ). These components are estimated by the mobile terminal by means of pilot symbols transmitted by the base station in the different directions θ_(k). The components are estimated periodically by the mobile terminal and transmitted, via the uplink channel, to the base station. At each new estimation, the integration calculation is reinitialised with the new components transmitted.

[0103] The vector {overscore (C_(d))} having been estimated by (26), the optimum weighting vector on transmission is then obtained by means of equation (21) in which the noise power matrix D_(d) is estimated by means of (23) and (24).

[0104]FIG. 4 illustrates an example of a device according to one embodiment of the invention. For reasons of simplicity, the processing of only one communication with a mobile terminal has been shown. The device, installed at the base station, comprises an array of antennae (40 ₀),(40 ₁), . . . (40 _(N−1)) coupled by means of duplexers to a first weighting module (41 ₁) on reception, weighting the signals received by the different antennae by means of a first weighting vector, {overscore (b_(u))}, and a second weighting module (41 ₂) on transmission, weighting a signal to be transmitted by a second weighting vector, {overscore (b_(d))}. When the device manages several conmmunications with a plurality of mobile terminals, other weighting modules identical to the modules (41 ₁), (41 ₂) must be provided in parallel with the latter. The signals received by the different antennae are demultiplexed according to the different users (that is to say the different mobile terminals). If the telecommunication system is a system of the DS-CDMA type, the demultiplexing will be carried out for example by means of a battery of filters adapted to the signatures of the different users, denoted (49). The N outputs of adapted filters are directed to a channel former (45), forming M beams in the directions θ_(k), k=0, . . . , M−1. The M signals are transmitted to a module (44 ₁) for evaluating the noise power matrix D. and to a module (43) for estimating the (uplink) channel {overscore (C_(u))}. Advantageously, the vector {overscore (C_(u))} is estimated by means of pilot symbols transmitted by the mobile terminal. Advantageously, the noise power matrix D_(u) is estimated during periods of silence of the mobile terminal. The matrix D_(u) and the vector {overscore (C_(u))} are supplied to a module (42 ₁) which calculates the vector {overscore (b_(u))} according to equation (16) and transmits it to the weighting module (41 ₁). The vector {overscore (C_(u))} is then differentiated by a differentiating filter (46). The resulting vector {overscore (ΔC)}_(u) is multiplied at (47) by the matrix M in order to give the vector {overscore (Δ)}_(d). This vector then servers for integrating the vector {overscore (C_(d))} in the integrator (48). This integrator is regularly initialised by the estimations of the downlink channel transmitted by the mobile terminal, denoted {overscore (C_(d))}(0). The vector {overscore (C_(d))} is transmitted to the matrix-type calculation module (42 ₂). This module also receives from (44 ₂) the noise power matrix D_(d). This matrix is evaluated in (44 ₂) by means of equation (23). To do this, the module (44 ₂) receives an estimation of the coupling coefficient, γ_(d) or Γ according to circumstances, of the interference powers σ_(dk) ² in the directions θ_(k) as well as the total power I_(d). The values σ_(dk) ² are advantageously calculated from equation (24) using the values of the transmission signals, S_(d)(j), j≠j₀, intended for the mobile terminals other than the one in question (j₀) and the gain vectors, {overscore (G)}_(d)(j), j≠j₀, which are associated with them. The module (42 ₂) effects the calculation of the vector {overscore (b_(d))} in accordance with equation (21) and transmits it to the weighting module (41 ₂).

[0105] Although the device described above has been depicted schematically in the form of functional modules, it goes without saying, however, that the various functions executed may be executed by virtue of a processor programmed for this purpose or by a plurality of dedicated processors. 

1. Method of obtaining a transmission gain function by means of an array of antennae, a signal to be transmitted by the array being weighted by a vector ({overscore (b)}_(d)) of N complex coefficients, referred to as the transmission weighting vector, N being the number of antennae in the array, the array transmitting to a telecommunication terminal over a transmission channel, referred to as the downlink channel, a downlink transmission signal (S_(d)) and the said terminal transmitting to the said array over a transmission channel, referred to as the uplink channel, an uplink transmission signal (S_(u)), the said downlink channel being disturbed by an isotropic noise (N′) and/or a directional noise, referred to as the downlink interference (I_(d)), characterised in that the said transmission weighting vector ({overscore (b)}_(d)) is determined by means of a matrix product from a noise power matrix (D_(d)) which is a function of the power of the said isotropic noise and/or of the power of the said directional noise and a vector ({overscore (C_(d))}), referred to as the downlink channel vector, representing an angular sampling of the transfer function of the downlink channel in M directions k, k=0, . . . , M−1, belonging to the angular range covered by the array.
 2. Method of obtaining a transmission gain function according to claim 1, characterised in that the said downlink channel vector ({overscore (C_(d))}) is obtained from variations in the transfer function of the uplink channel.
 3. Method of obtaining the transmission gain function according to claim 2, characterised in that the said downlink channel vector ({overscore (C_(d))}) is obtained from variations ({overscore (ΔC_(u))}) in a vector ({overscore (C_(u))}), referred to as the uplink channel vector, representing an angular sampling of the transfer function of the uplink channel in the said M directions.
 4. Method of obtaining a transmission gain function according to claim 3, characterised in that the variations ({overscore (ΔC_(d))}) in the downlink channel vector are obtained from variations ({overscore (ΔC_(u))}) in the uplink channel.
 5. Method of obtaining a transmission gain function according to claim 3 or 4, characterised in that the variations Δc_(dk) in the components c_(dk) of the downlink channel vector ({overscore (C_(d))}) are obtained by means of the variations Δc_(dk) in the components c_(uk) of the uplink vector by: Δc_(dk)/c_(dk)=f_(d)/f_(u).Δc_(uk)/c_(uk) where f_(u) is the frequency used on the said uplink channel and f_(d) is the frequency used on the said downlink channel.
 6. Method of obtaining a transmission gain function according to claim 4 or 5, characterised in that the said downlink channel vector ({overscore (C_(d))}) is obtained by integrating the said variations ({overscore (ΔC_(d))}) in the said downlink channel vector and an initial value ({overscore (C_(d))}(0)) transmitted by the said terminal.
 7. Method of obtaining a transmission gain function according to one of the preceding claims, characterised in that the noise matrix is a diagonal matrix of size M×M and of components {square root}{square root over (σ_(dk) ²+γ_(d)N′₀/I_(d))} where σ_(dk) ² is the power of the downlink interference in the direction k, N′₀ is the power of the isotropic noise, γ_(d)=1/∥{overscore (C)}_(d)∥² and I_(d) is the total power of the downlink interference.
 8. Method of obtaining a transmission gain function according to one of claims 1 to 6, characterised in that, the array transmitting over a plurality of downlink channels a plurality of transmission signals to a plurality of telecommunication terminals and receiving from them a plurality of transmission signals transmitted over a plurality of uplink channels, each downlink channel j relating to a terminal j of the said plurality being associated with a transmission weighting vector {overscore (b)}_(d)(j), the second noise matrix relating to the downlink channel j is a diagonal matrix of size M×M and of components {square root}{square root over (σ_(dk) ²(j)+γ_(d)(j).N′₀/I_(d)(j))} where σ_(dk) ²(j) is the power of the downlink interference for the downlink channel j in the direction of k, γ_(d)(j) is a coefficient characterising the power transfer over the downlink channel j, N′₀ is the power of the second isotropic noise, and I_(d) is the total power of the downlink interference.
 9. Method of obtaining a transmission gain function according to claim 8, characterised in that the coefficient γ_(d) (j) is transmitted to the array by the terminal j on the associated uplink channel.
 10. Method of obtaining a transmission gain function according to claim 8, characterised in that the coefficient γ_(d) (j) is estimated by the base station from a coefficient (r) characterising the power transfer in the uplink direction.
 11. Method of obtaining a transmission gain function according to one of claims 8 to 10, characterised in that, for a given downlink channel j, the downlink interference power in the direction k, σ_(dk) ² (j), is estimated according to the power of the transmitted signals (S_(d)(j′)) on the downlink channels j′ distinct from j, a coefficient β_(d)(j) characterising the orthogonality of the downlink channel j, the components g_(dk)(j′)) of the gain vectors ({overscore (G)}_(d)(j⁻¹)) relating to the said distinct downlink channels j′ the gain vectors consisting of an angular sampling in the said M directions of the transmission gain functions obtained for the said distinct downlink channels j′.
 12. Method of obtaining a transmission gain function according to claim 11, characterised in that the said coefficient β_(d)(j) is estimated from a coefficient characterising the orthogonality of the uplink channel j.
 13. Transmission device for a base station in a mobile telecommunication system, comprising an array (40 ₀,40 ₁, . . . ,40 _(N−1)) of N antennae, weighting means (41 ₂) for weighting the signal to be transmitted (S_(d)) by the said array by means of a transmission weighting vector ({overscore (b)}_(d)) of N complex coefficients, characterised in that it comprises means (42 ₂,44 ₂,46,47,48) adapted to implement the method of obtaining the transmission gain function according to one of the preceding claims, the said adapted means supplying to the said weighting means (41 ₂) the said transmission weighting vector ({overscore (b)}_(d)). 